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We study spin dynamics for two electrons confined to a double quantum dot under the influence 
of an oscillating exchange interaction. This leads to driven Rabi oscillations between the |t|)-state 
and the m)-state of the two-electron system. The width of the Rabi resonance is proportional 
to the amplitude of the oscillating exchange. A measurement of the Rabi resonance allows one to 
narrow the distribution of nuclear spin states and thereby to prolong the spin decoherence time. 
Further, we study decoherence of the two-electron states due to the hyperflne interaction and give 
requirements on the parameters of the system in order to initialize in the |tJ.)-state and to perform 
a VSWAP operation with unit fidelity. 



CD 



I 

O 

o 



> 



o 



I 

O 

o 



X 



PACS numbers: 73.21. La,76.20.+q,76.30.-v,85.35.Be 



I. INTRODUCTION 

One of the important proposals for quantum infor- 
mation processing in solid-state systems is the spin- 
qubit proposal for quantum computing with electron 
spins in quantum dotsi. Much effort has been put 
into the realization of this proposal leading to exciting 
theoretical^ and experimental achievements li^'^'SiLSiS 
Still many challenges remain such as decoherence and 
the implementation of single-qubit gates. 

A major obstacle to quantum computation with the 
quantum-dot spin qubit is decoherence due to the cou- 
pling of the qubit to its environment. The hyper- 
fine interaction between the electron spin and the nu- 
clear spins present in all III-V semiconductorsiS leads 
to the strongest decoherence effect^i"i^^i^'^i^'^i^^i^^i^'^i^^. 
ExperimcntsSiSiiSiSS have yielded values for the free- 
induction spin dephasing time T| that are consistent 
with Tl - VN/A ~ lOnsili^ for TV = 10^ and 
A = 90^eV in GaAs, where TV is the number of nu- 
clei within one quantum dot Bohr radius and A char- 
acterizes the hyperfine coupling strengthi^ This is to be 
contrasted to potential spin-echo envelope decay, which 
may be much largeri22iSi^ With a two-qubit switching 
time of Ts ^ SOpsii this only allows ^10^ gate opera- 
tions within T2, which falls short (by a factor of 10 to 
10^) of current requirements for efficient quantum error 
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There are several ways to overcome the problem of 
hyperfine-induced decoherence, of which measurement 
and thus projection of the nuclear spin state seems to 
be the most promising onaii. Other methods include 
polarizationii4iSiiL2& of the nuclear spins and spin echo 
techniques i2iiL2^ However, in order to extend the de- 
cay time by an order of magnitude through polariza- 
tion of the nuclear spins, a polarization of above 99% 
is required^ but the best result so far reached is only 
~60% in quantum dotsi^i^ With spin-echo techniques, 
gate operations still must be performed within the single- 



spin free-induction decay time, which requires faster gate 
operations. A projective measurement of the nuclear spin 
state leads to an extension of the free-induction decay 
time for the spin. This extension is only limited by the 
ability to do a strong measurement since the longitudi- 
nal nuclear spin in a quantum dot is expected to survive 
up to the spin diffusion time, which is on the order of 
seconds for nuclear spins surrounding donors in GaAs^ 

The implementation of quantum computation schemes 
requires coherent control of the qubits. Rabi oscillations 
between the two qubit states arc an important signature 
of coherence and thus observation of controlled Rabi os- 
cillations is an important intermediate step in the exper- 
imental implementation of quantum information proces- 
sors. Despite recent experimental achievements;^ there 
has still been no experimental observation of driven Rabi 
oscillations for a system of two quantum-dot spin qubits. 
What has been observed is electron spin resonance via g- 
tensor modulation in a bulk semiconductor 1^ 

In the quantum-dot spin qubit proposal, two-qubit 
gates are realized through tuning of the exchange cou- 
pling J between the two spinsiii The splitting between 
singlet and triplet states of the two-electron system is 
given by the exchange coupling J and in devices such 
as those in Refs. and H, J can be controlled through 
gate voltages. Petta e t al^ have recently managed to 
implement the V S WAP-gate in their setup. However, in 
order to implement single-qubit gates, control over local 
magnetic fields or g-factors is required^ 

As we will show in Sec^l an oscillating exchange J(i) 
induces Rabi oscillations between the states |ti) and ||t) 
of two electron spins (one electron in each dot) . The am- 
plitude of these oscillations is resonant on the splitting 
between ||J,) and | jf) and the width of this resonance is 
proportional to the amplitude j of the oscillating com- 
ponent of J{t) = Jq + jcoslujt), where lu is the driving 
frequency. Since the splitting depends on the state of 
the nuclear system, a measurement of the resonance is 
also a measurement of the state of the nuclear spins and 
thus provides a way to narrow the quantum distribu- 



2 



tion of the nuclear spin states. This narrowing of the 
spin state is one possible solution to suppress hyperfine- 
induced decoherence in quantum-dot spin qubita^. It 
has been proposed to measure the nuclear spin polar- 
ization using a phase estimation methodj^S In the ideal 
case, phase estimation yields one bit of information about 
the nuclear-spin system for each perfectly measured elec- 
tron. Optical methods have also been proposedi^ The 
all-electrical method we present here can be applied with 
current technology. 

The rest of this paper is organized as follows. In Sec. |nl 
we show that an oscillating exchange leads to driven Rabi 
oscillations and calculate the resonance linewidth. In Sec. 
Illll we propose a method to narrow the distribution of the 
nuclear spin states, in Sec. II VI we consider decoherence 
induced through the hyperfine interaction for a static ex- 
change coupling J. Wc use these results in Sec|V]to ana- 
lyze under which conditions we reach unit fidelity for the 
initialization to the state jtl) and a v&WAP operationi 
Sec. IVII contains a summary of our results. 

II. OSCILLATING EXCHANGE AND ESR 

In this section we show that under suitable conditions 
an oscillating exchange interaction may be used to induce 
Rabi oscillations in a system of two electrons confined to 
a double quantum dot like those in Refs. I^QIIS 

We denote by hi = (/if , /if , /if ), i = 1,2, the col- 
lective quantum nuclear spin operator, the "Overhauser 
operator" , in dot one and two, respectively, and write 
Sh^ = \{hl — /12). The collective quantum nuclear spin 
operator is defined as = X^fe^fc-'-fe' where I^. is the 
nuclear spin operator for a nucleus of total spin / at lat- 
tice site k, and the hyperfine coupling constants are given 
by A], = wv4|'(/'o(rfc)P, where v is the volume of a unit cell 
containing one nuclear spin, A characterizes the hyper- 
fine coupling strength, and ipl{rk) is the single-particle 
envelope wavefunction of the electron evaluated at site 
k. Further, {0)rms = (V'/l l^i)^^^ is the root -mean- 
square expectation value of the operator O with respect 
to the nuclear spin state \ipi)- We assume that the Zee- 
man splitting Cz = gfJ-sB induced by a uniform applied 
magnetic field B = (0,0, 5), B > 0, is much larger than 
((5h)i.nis and (hi)i.,„s. Under these conditions the relevant 
spin Hamiltonian becomes block diagonal with blocks la- 
beled by the total electron spin projection along the mag- 
netic field 5*^ . In the subspace of 5^ = the Hamiltonian 
can be written as {h — l)i£ 

Ho^^{l+T^)+Sh^T^ +db'T''. (1) 

Here, J is the Heisenberg exchange coupling between 
electron spins on the two dots and Sb^ the inhomogene- 
ity of an externally applied classical static magnetic field 
which we add in addition to the treatment in Ref. El 
Further, r = (T^,T^,r^) is the vector of Pauli matri- 
ces in the basis of S"^ = singlet \S) and triplet |To) 



(IS") ^ |t^ = -1) , |To) = +1)). It has been pro- 

posed to use two pseudo-spin states such as \S) and \Tq) 
as a logical qubiti^ 

We assume a time-dependent exchange of the form 

J = J{t) = Jo + i cos(wi). (2) 

The operator (5/i^ commutes with the Hamiltonian at all 
times. Thus, if the nuclear-spin system is in an eigen- 
statc |ri) of Sh^ with Sh^ \n) = Shf^ \n), we have H \tp) = 
Hn iV'e) ^ where in i/„ the operator 6h^ has been 
replaced by dh^ and \^e) is the electron spin part of the 
wave function. In order to bring i/„ to a form that is very 
similar to the standard ESR (electron spin resonance) 
Hamiltonian^ [Hesr = — ^CzCz — \^^xCOs{ujt)ax) we 
perform a unitary transformation Ui = exp(— i-|T^) 
which is just a rotation about the y-axis in a Bloch- 
sphere picture. Also introducing ri„ = 2{5h^ + 5h^), the 
above Hamiltonian becomes 

= UiHnUl ^^T-' + l cos{ujt)T- - hl^r'. (3) 

The Pauli matrices are now given in the new basis of 
lit) = = 1) = 1+) and ITi) = |r- = -1) = |-). For 
Jo = this is just the standard ESR Hamiltonian. We 
have evaluated pseudo-spin dynamics under this Hamil- 
tonian in a rotating wave approximation close to res- 
onance for j ^ r2„. When we treat the Jo-term as 
a perturbation and calculate the transition probability 
between unperturbed eigenstates of the Hamiltonian we 
find that it is proportional to Jq/17;:^ and we may thus ne- 
glect this term close to resonance and if Jq <C iln- Hence, 
we are left with the standard ESR Hamiltonian which 
leads to Rabi oscillations. Initializing the two-electron 
system in the state |J,|) = |+) (which can be done as 
proposed in Sec. 0) we obtain for the expectation value 
of r^(<): 

(r-(i))„ - {n\®{+\T'{t)\+)®\n) 

J = 2{(si~^jf^{m\ (5) 

j < r2„, Jo < r2„, - w| < r2„. (6) 

For 10 = f2„ the system undergoes coherent Rabi oscilla- 
tions between the states |-|-) and |— ) with a frequency of 
j. Averaged over time, the expectation value of is 

((r^)„) = hm 1 C {T\t))^dt = 

T-»oo T Jo ("n - + (.7/2)^ 

(7) 

In order to measure the time-averaged value ((t")„) the 
measurement time must be much larger than the period 
of Rabi oscillations (~ 1/j on resonance). 1 — ((t^),i) 
has a Lorentzian lineshape with a full width at half maxi- 
mum (FWHM) of J. Most importantly, the resonance fre- 
quency depends on the nuclear-spin eigenstate through 
51„ — 2((5/i^ -I- 56^) and thus a measurement of the reso- 
nance will determine bh/L . 
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A. Superposition of nuclear— spin eigenstates 

Before a measurement on the nuclear-spin system is 
performed, there is no reason for the nuclear-spin sys- 
tem to be in an eigenstate of 6h', but it is most likely 
in some generic superposition of these eigenstates. Thus, 
we now investigate how the resonance changes if we con- 
sider the nuclear-spin system to be in a superposition of 
eigenstates of the collective nuclear spin operator (5/i^ . 

At t = we fix the electron system in the state |J,t) = 
|-|-) while the nuclear-spin system is in an arbitrary state; 
p{0) = pe(0) ® /9/(0) with 



Pe(0) 

P/(0) 



a„ n 



(8) 
(9) 



where the a*j satisfy the normalization condition 
E„KP = 1 and Y.^P^ = 1- Here, pi{n) = Y.^P^K\^ 
are the diagonal elements of the nuclear-spin density op- 
erator. The Hamiltonian Hq commutes with 5h^ and 
thus we find 



(10) 



which defines the overbar. 

We assume that for a large number of nuclear spins 
^ 1 which are in a superposition of (5/i^-eigenstatcs 
pi{n) describes a continuous Gaussian distribu- 
tion of 5h^ values, with mean 5h^ and variance = 



{5h^ — Sh^^ . In the limit of large N the approach to 
a Gaussian distribution for a sufficiently randomized nu- 
clear system is guaranteed by the central limit theorem^ 
We perform the continuum limit according to 

^Pi{n)f{n) j dxpi.-x,a{x)f{x), (11) 
Pi::xm{x) = 



1 I ix — x)^ , , „, 

exp I I , (12) 
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where x = Shf^, x 



Sh^ and cr^ = x'^ — x'^. The only 
effect of 6b^ is to shift the mean value of the Ovcrhauser 
field inhomogencity to = x + 5b^ , whereas the width 
is left unchanged: = a. According to this description 
we obtain 



(r^(t)) 
g{x,t) 



dxpi-xo.aoix){f{x)+g{x,t)), (13) 
(2x-a;)2 



(2x-a;)2 + (j72)2' 
(j72)2 cos (2^{2x~lo)^ + {j/2)H) 



{2x - Ljf + {j/2Y 



(14) 
(15) 



The second term fEa. HlSf) ') vanishes when it is averaged 
over time and wc find 



J 



1 



2<TnV27r 



dx exp 



(x - 2.To)^ 



(.?72)^ 



072)^ 



(16) 
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This integral (a convolution of a Lorentzian and Gaus- 
sian) is the well-known Voigt functionjS and the result- 
ing lineshape is the so-called "Voigt profile" . The Voigt 
function may be expressed as (a) = j + 41to — 2iLo) 



= 1 - 



4crn 



exp 



32crg 



erfc 



(4V2, 



(17) 



where erfc(z) is the complementary error function. In the 
regime where ctq j we may approximate the Lorentzian 
in the convolution fEa. H16|) )bv its value at x = 2xq and 
obtain 



{2x(i~ujf 



(2xo-c.)2 + (j72)2' 



(18) 



In this case the resulting resonance has the same FWHM 
as the Lorentzian, viz. j. On the other hand, if (Tq 3> j, 



we may approximate the Gaussian with its value a,t x — u 
and thus obtain 



J 
4(70 



2 exp 



(2x0 - 



CTO > j. 



(19) 

In this regime the width is twice the width (Tq of the 
Gaussian distribution of the nuclear spin states. In order 
to make a statement about the width of the Voigt profile 
in general we look at the peak-to-peak separation Ay of 
the first derivative of the Voigt profile. For a Gaussian 
with a standard deviation of 2cro wc find Ag = 4cro for 
the peak-to-peak separation of the derivative and for 
a Lorentzian with FWHM of j we have Ai, = j /y/2>. A 
Fade approximant for Ay in terms of A^ and Ag yielda^i 



Av 



Ag; + oiAgAl + asAl 
Ag + a2AL 



(20) 
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where ai = 0.9085, 02 = 0.4621. This approxima- 
tion is accurate to better than O.OlAy for all values 
of Al, a similar formula may also be given for 

the half width at half maximum (HWHM) of the Voigt 
profile 



III. STATE NARROWING 

The general idea behind state narrowing is that the 
evolution of the two-electron system is dependent on the 
nuclear spin state and thus knowing the evolution of the 
two-electron system determines the nuclear spin state. 
Thus, in this section we describe how the Gaussian super- 
position pi-ao,xo{x) of collective nuclear spin eigenstates 
\n) can be narrowed through a sequence of measurements 
performed on a double quantum dot on a time scale much 
less than the timescale of variation of 6h^ and for j < cro- 
We first give a general description of how a complete mea- 
surement of the lineshape of the Rabi resonance narrows 
the Gaussian superposition. Such a complete measure- 
ment of the lineshape consists of many single measure- 
ments of the operator r^. In Sec. IIII AI we present a 
detailed analysis of such a complete measurement and in 
Sec. liiilji we discuss different measurement schemes. 

The operator Sh^ was defined in Sec. HTland it describes 
the difference in the z-components of total nuclear field in 
each of the two dots. The total nuclear field is the result 
of ~ 10^ single nuclear spins and thus the eigenvalues 
of Sh^ will be highly degenerate. In the limit of large TV 
the spectrum of Sh'^ is quasi-continuous and the proba- 
bility density of eigenvalues of Sh^ is given by a Gaussian 
distribution, as described in Sec. Ill AI For such a Gaus- 
sian superposition of nuclear spin eigenstates, the line- 
shape of the Rabi resonance is given by a Voigt profile, 
as described in Sec. Ill AI This Voigt profile can be seen 
as a superposition of Lorentzian lineshapes, where each 
Lorentzian results from a nuclear spin eigenvalue Shf^ and 
is centered around il„ = 2{6hf^ + 6b^). In the Voigt pro- 
file, these Lorentzian lineshapes are weighted according 
to the amplitude of the corresponding eigenvalue Shf^ in 
the Gaussian-distributed superposition. Through a per- 
fect complete measurement of the Rabi-resonance line- 
shape, the superposition of Lorentzian lineshapes col- 
lapses and we are left with one single Lorentzian (see fig- 
ure^. This Lorentzian corresponds to one single eigen- 
value of dh^ and thus the Gaussian distribution has been 
narrowed to zero width; the nuclear-spin system is in a 
state with fixed eigenvalue Sh^. 

In principle, we would need to do infinitely many single 
measurements in order to completely measure the line- 
shape of the Rabi resonance with perfect accuracy, since 
each point on this resonance curve is a (time-averaged) 
expectation value of the quantum mechanical operator 

. Still, we may perform a finite number M of single 
measurements (see Sec. IIII A|) for each of a set of driving 
frequencies u; and thus obtain the series of expectation 
values for different uj up to some error. This error de- 
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FIG. 1: a) This figure illustrates the projection ob- 
tained through an ideal complete measurement of the Rabi- 
resonance lineshape. All the different Lorentzian resonances 
corresponding to different nuclear spin eigenstates add up to a 
Gaussian lineshape. b) Through a perfect complete measure- 
ment of the lineshape of the Rabi resonance, which involves 
many single measurements of r^, the superposition collapses 
and we are left with one single Lorentzian centered around 
2x'q = Qn, which in general is different from 2x0. 



pends on AI . There will then in general be more than 
one Lorentzian which can be fit (within error) to these 
expectation values and thus we would not narrow to zero 
width. We would still have a distribution of nuclear spin 
eigenstates, but one with smaller width than before the 
measurements. 

For such a narrowing through measurement to be suc- 
cessful, the amplitude j of the oscillating exchange J{t) 
which determines the width of the Lorentzian lineshapes 
should be smaller than the width do of the Gaussian dis- 
tribution. Otherwise, the Rabi resonance would be dom- 
inated by the Lorentzian (see Eq.l^SJ) and the method 
would not result in narrowing of the nuclear-spin distri- 
bution. The general requirements on the system parame- 
ters to narrow the distribution of nuclear spin eigenvalues 
arc 



j,Jo,cro < xo; 



3 ^ f^o- 



(21) 
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We note that, unlike in standard ESR, power absorption 
is not measured here, but instead the expectation vahie 
of the pseudo-spin , for instance via a quantum point 
contact (QPC) nearby one quantum dot (for a detailed 
description of the measurement process via such a QPC 
we refer the interested reader to Ref. Isfil) . To determine 
the expectation value of the pseudo-spin many single 
measurements of the pseudo-spin are necessary and we 
thus proceed to give a detailed description of the state 
narrowing by considering the effect of these single mea- 
surements on the nuclear spin state. 



we thus assume^ that measurements are performed with 
low time resolution At ^ Further, we perform the 
continuum limit and obtain for the probability density of 
eigenvalues, i.e., the diagonal part of the density operator 
in the continuum limit (with x = 5h^ + Sb^ and pi{x) = 
Pi;xo,ao{x), see Eq.lO): 

p^^'+^-)(x) = p,(x)(l-i„(x))-^, (25) 



(26) 



A. Description of state narrowing by consecutive 
pseudo— spin measurements 

In this subsection we describe in detail how a single 
measurement of the pseudo-spin of the two-electron 
system affects the nuclear-spin system. Further, we 
give a general formula for the diagonal elements of the 
nuclear-spin-system density operator in the continuum 
limit after M measurements. The sequence of M mea- 
surements is referred to as a "complete measurement" . 

At t = the two-electron system is initialized to the 
state |-|-) = \\X) and we assume that the electron and the 
nuclear system are initially factorized. Thus, the total 
system at t = is described generally by the following 
density operator 



p(0) = pe(0)®p/(0) 



(22) 



with nuclear-spin state = X^n '^n 1'^)- The diagonal 
elements of the nuclear-spin density operator at t = are 
given by pi(n) = pi(n,0) = X^j-PiknP ^^"^ ™ the con- 
tinuum limit we obtain the probability density pi;x,a{x) 
for the eigenvalues (5/i^ = a; as given in Ea. H12|l . At time 
tm a measurement of the two-electron system (at driving 
frequency uj, where to is defined in Eq.(|2Jl) is performed 
with two possible outcomes |-|-) and |— ). The diagonal 
elements of the nuclear-spin density operator after the 
measurement are given by (see Appendix IB|| 

(1,±), , N P/('^70)l, , 2, , , , 

P\ '{n,trn)^ p^^^ ^ -(1±(t {trn))n), (23) 

where (T^(i))„ is given by Eq.Q) and the probabilities 
P^{tm) to measure |±) are 



P±(t,„) = i (1 ± (r^(i„0)„)p,|< 



(24) 



In the case where a measurement is performed with a 
low time resolution^ At, i.e., if At ^ the density 
operator after the measurement is the time average over 
the time interval At and the cosine term in {T'^{tm))n 
averages out (note that in the case of a measurement 
with low time resolution, tm is arbitrary, as long as At is 
chosen to be large enough) . For the rest of this subsection 



where the probabilities for measuring |-|-) or |— ) are given 

by 



P+ 



dxpi{x){l~L^{x)), (27) 
dxpi{x)L^{x), (28) 



with 



L^{x) 



(j74)= 



2 (x- 



w ^2 
2 ) 



074)2 



(29) 



After the first measurement, the two-electron system 
is reinitiahzed to the state |-|-) if necessary and a sec- 
ond measurement is performed. Since the initial den- 
sity matrix factors out in the above results, it is clear 
how to generalize Eas. (|25|l and H26|) to the case where 
A/ consecutive measurements (without randomization of 
the nuclear-spin system in between measurements) are 
performed: every time |-|-) is measured, the diagonal el- 
ements pi{x) of the nuclear density matrix is multiplied 
by 1 — L^ix) and every time |— ) is measured, pi{x) is 
multiplied by L^{x). Thus, we obtain the diagonal el- 
ements of the nuclear density matrix after 
Af measurements, of which a~ times the measurement 
outcome was |— ) (and {M — Q!~)-times |-|-)): 



(M,a- 
Pl 



Pl{x) 



-Wuj{M, a~]x). 



(30) 



Here, Wuj{M,a ;x) and the normalization factor 
Quj{M, a~) are given by 

WUM,a-;x) = L^(a;)"" (1 - L^(a;))^^-"^ (31) 
Q^{M,a-) = / dxpi{x)W^{M,a-;x). (32) 



The normalization factor Q^{M,a ) is related to 
through P- = gc^(l,l), P+ = QJil,0). In the case 
where measurements are performed at rrif different fre- 
quencies, Ea. l|Sn|) generahzes to 



(33) 
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The probability density ^^^'^'^^ (a;) after M 

measurements performed at different driving frequen- 
cies depends on the frequencies {wj} = {wi, . . . , Wm/}, 
the number of measurements at each frequency {Mi\ = 
{Ml, . . . , Mm;}, and the number of times |— ) was mea- 
sured at each frequency {a~} = {oti-, . . . , ot^ }. Ea. H33|) 
gives the distribution of nuclear spin eigenvalues for any 
sequence of M measurements, i.e., without randomiza- 
tion of the nuclear-spin system in between measure- 
ments. 



B. Measurement schemes 

In this subsection we describe different measurement 
schemes. One main characteristic of the schemes is 
whether we have unconditional evolution of the nuclear- 
spin density matrix between measurements (one waits 
for the nuclear-spin system to rerandomize between sub- 
sequent measurements), or whether we have conditional 
evolution, i.e., the nuclear-spin system is assumed to be 
static between measurements. 



1. Unconditional scheme 

The simplest scheme is to measure only once at one 
single driving frequency u). If the outcome is |— ), the 
nuclear-spin distribution after the m easurement is given 
by Eq.(j2ni); the FWHM (2croV2 In 2 « 2cro) of the initial 
distribution will have been narrowed by a factor « j/4(To 
(the nuclear-spin distribution will approximately be a 
Lorcntzian with FWHM of j/2). For j <C o-q and 
Lo = 2a;o, the probability P~ to measure |— ) in the first 
measurement is P^^2xa ~ i/6fo (the exact formula is 
given in Eq.jSSJ). If the measurement outcome is |— ), 
we stop measuring. Otherwise, we wait for the system to 
rerandomize (in contrast to the conditional schemes) and 
perform another measurement. This is repeated until |— ) 
is measured for the first time. On average one needs to 
perform Ad' w Gctq / j measurements in order to narrow by 
a factor of w j/icro (we write M' because this number of 
measurements should not be confused with the number of 
measurements M used above in the case of measurements 
performed without rerandomization in between). If the 
driving frequency lu is far from the center xq of the ini- 
tial Gaussian distribution, the number of required mea- 
surements increases by a factor of exp((a;o — a;/2)^/2crQ). 
This always leads to a narrowed distribution which is 
centered around lo/2. Thus, with this scheme it is possi- 
ble to choose the center of the nuclear-spin distribution 
after the measurement. This unconditional measurement 
scheme is the one which should be easiest to implement 
in an experiment since one only needs to measure once at 
one single frequency. However, if measurements at sev- 
eral diiferent frequencies can be performed, a systematic 
narrowing of the distribution can be implemented as we 
show next. 



2. Adaptive conditional scheme 

The probability of measuring |— ) in a measurement is 
determined by the overlap of the Lorcntzian L^{x) and 

the probability density of eigenvalues P/^^'" (for 
the first measurement this probability is P~ , which is 
given in Ea. (|28|l 'l. Then, if we have the outcome |— ) 

for a measurement at driving frequency w, pf^'" 

as a function of x becomes peaked around uj/2 ( since 

L^{x) is centered around x = w/2), the overlap of the 

Lorcntzian L^^{x) and p\^^'°' '"^(a;) increases and there- 
fore the probability to measure |— ) in a subsequent mea- 
surement also grows. If, on the other hand, we have out- 
come |-f), the term l~Li^{x) causes a dip in p\^^'" 
at X = ij-'/2, the overlap of the Lorcntzian L^{x) and 

p\^^'" decreases and thus the probability to mea- 

sure |— ) in a subsequent measurement with the same 
driving frequency uj also decreases. Since it is the mea- 
surement outcome |— ) that primarily leads to narrowing, 
the measurement scheme should maximize the probabil- 
ity to measure |— ). This can be achieved by changing 
the driving frequency to always in such a way that be- 
fore each measurement L^{x) and the nuclear-spin dis- 
tribution pj*^'" ''^\x) have their maximum at the same 
x, i.e., set a;/2 = a;max, where a;,„ax is the x for which 
^^M,Q '^^ [x) has a maximum. Thanks to the adaptive 
driving frequency lo, the probability P^J to measure |— ) 
is « j /6cFo in each measurement until |— ) is measured for 
the first time. Without adapting, i.e., when measuring 
always at the same driving frequency w, P~ decreases, 
as explained above (as long as we do not measure |— )). 
After measuring |— ) for the first time, the probability P~ 
to measure |— ) increases. Every time the measurement 

outcome is |— ), the distribution pj*^'" ''^\x) is multi- 
plied by L^{x) and becomes narrower (since L^{x)°' 
has a FWHM of {j/2)V2^/'^' - 1). However, the mea- 
surement outcome |+), for which pj*^'" '"\x) is multi- 
plied by 1 — L^{x), is still more likely and leads to a 
small widening of the distribution. Our simulations of 
this measurement scheme do, however, show that after 
|— ) has been measured several times, the nuclear spin 
distribution is narrowed by more than a factor j/4cro- 

This adaptive scheme was first proposed in an optical 
setup by Stepanenko et al. in Ref. [S^. This scheme re- 
quires that Xniax can be calculated (or read from a table) 
between subsequent measurements and that the driving 
frequency uj can be tuned with a precision that is bet- 
ter than the width of the nuclear-spin distribution be- 
fore each measurement. For this adaptive scheme (and 
other conditional schemes) to work, it is important that 
the nuclear-spin system does not randomize during the 
course of the complete measurement, i.e., the complete 
measurement must be carried out within a time that is 
shorter than the time scale for nuclear spin dynamics. We 
thus assume that the nuclcar-spin system (viz. (5/i^) has 
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no internal dynamics between the single measurements 
of (t) , but only changes due to the measurements per- 
formed on the two-electron system, i.e., due to single 
measurements of t^(<). We expect Sh^ to vary on the 
time scale of nuclear spin diffusion out of the dot, which 
is on the order of seconds for nuclear spins surrounding 
donor impurities in GaAs^ However, there may be other 
sources of nuclear spin dynamics (see also Appendix 
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FIG. 2: In this figure we show a typical— se- 
quence of the rescaled probability density of eigenvalues 



for 



7r(j;) = pj {x)/ma.xlpj ( 

the adaptive conditional scheme. Here, p^''*^»b{aj ^'^'^^^^(^x) 
is given in Eg.ll^. We have x = 5h^ + 5V , j/uo = 1/10 
and in a)-c) the initial Gaussian distribution (with FWHM 
2(Jo V2 In 2 ~ 2aa) is plotted for reference, a) Up to M = 50 
measurements the outcome is never |— ) and thus each mea- 
surement "burns a hole" into the distribution where it previ- 
ously had its maximum, b) In the SI''* measurement the out- 
come is |— ) which leads to a narrowed distribution of nuclear 
spin eigenvalues (peak centered at « 0.5) with a FWHM that 
is reduced by a factor ~ j/4o-o- c) Adapting the driving fre- 
quency uj to this peak, i.e., setting uj/2 = a::max in subsequent 
measurements, leads to further narrowing every time |— ) is 
measured. In this example the final FWHM is « cto/IOO, i.e., 
the distribution has been narrowed by a factor « j'/lOcro. d) 
The probability P~ to measure |— ) jumps up after the 51°' 
measurement and after |— ) is measured several more times, 
this probability saturates close to 1/2. 

In figure|21we show a typical^ sequence of nuclear spin 
distributions for the adaptive scheme with total number 
of measurements M = 100 and j /(Tq = 1/10. We see ( 
figure 121 (a)) that up to M = 50 the measurement out- 
come is never |— ) and thus each measurement "burns 
a hole" into the distribution where it previously had its 
maximum. In the 5P* measurement (figure|2Ib)) the out- 
come is |— ), which narrows the distribution by a factor of 
~ j/4:ao- Adapting the driving frequency lu to this peak, 
i.e., setting ui/2 ~ Xmax hi subsequent measurements, 
leads to further narrowing, i.e., to a total narrowing by 
more than a factor j/4tTo (figure l^^c)). In this example 
we have a~ = 22 after M — 100 measurements and the 



final FWHM is « ao/lOO, i.e., the distribution has been 
narrowed by a factor « j/lOtro- In figure [2Id) the prob- 
ability to measure |— ) before each measurement is 
shown. After the first time |— ) is measured, P~ jumps 
up and after several more times |— ) was measured, it 
saturates close to 1/2. P~ is a good signature of the dis- 
tribution's width. As the width of the distribution goes 
to zero, P^ approaches 1/2. This adaptive conditional 
scheme is more intricate than the unconditional scheme, 
but allows one to narrow by more than a factor j/4(To. 



3. Other conditional schemes 

Other possible measurement schemes involve measure- 
ments at several frequencies, as in the adaptive scheme. 
One may either choose a fixed number of frequencies 
within one or two do and measure several times at each 
frequency (without randomization between the measure- 
ments) or sweep the frequency, i.e., measure only once 
at each frequency but vary the frequency only in small 
steps. Based on numerical simulations of these schemes, 
we find that the typical number of measurements to nar- 
row by a factor of j /(Tq is greater than in the adaptive or 
the unconditional (single-frequency) schemes. 



4- Time-domain measurement scheme 

We note that when a complete measurement of one 
of the correlators discussed in Sec. IIVI is performed 
with perfect resolution in time and perfect accuracy, this 
would also determine the state of the nuclear spin system 
and thus narrow the distribution of nuclear spin states. 
This is because th e frequency o f the oscillating correla- 
tors is given by + ^{Sh^)^ and thus measuring the 
frequency of the correlator determines the eigenvalue Sh^ 
of the nuclear-spin system. However, it may be possible 
to perform a weak measurement of the decay of the corre- 
lators and thus also to see the prolongation of the decay 
after applying a narrowing scheme. To understand in de- 
tail the effect of measurements in the time domain, fur- 
ther study is required. Narrowing through measurement 
of the correlators is a time-domain measurement. In con- 
trast, the narrowing schemes we have proposed above are 
frequency-domain measurements. If the frequency res- 
olution is better than the time resolution, our method 
would most likely be more suitable. 



IV. CORRELATION FUNCTIONS IN THE 
= SUBSPACE 

In this section we investigate the Hamiltonian Hq of 
Eq. 1^ with static exchange couphng J. Using this 
Hamiltonian we wish to calculate correlation functions 
for several observables in the subspace of zero total spin 
in the z-direction. In our previous worki^ we calculated 
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the time evolution of a particular correlator involving the 
states \S) and \Tq). However, there are four additional 
independent correlators involving the x and y compo- 
nents of pseudo-spin which require a separate calcula- 
tion. Quite surprisingly, it will turn out that these corre- 
lators have different decay behavior in time. The correla- 
tors we calculate here show the decohcrence properties of 
the pseudo-spin states under the influence of the hyper- 
fine interaction. There may be additional sources of de- 
coherence which we do not consider here, such as orbital 
dephasing, corrections to the effective Hamiltonianfi^ the 
couphng of the QPC to the dot spins^^ etc. The results 
of this section will help to give requirements on the pa- 
rameters of the system in order to initialize in the state 
Iti) and to assess the fidelity of a vSWAP operation with 
static J (see Sec. Iv|l . 

Diagonalizing Hq gives the following eigenvalues and 
eigenvectors 



{nj2)\S} + Et\To 



«) \n) 



(34) 
(35) 



where again \n) is an eigenstate of the operator 5h^ with 
Sh" \n) = Shf^ \n). At t = we fix the electron system in 
an arbitrary superposition of \Tq) and 15) 



\Mt = 0)) = \A); \A) 



cos ^ 1 5)4-6'^- sin^lTo) 



(36) 

The nuclear-spin system is again in a general state (see 
Sec. Ill A|) . As will be shown in Sec. El it is possible, 
in principle, to initialize to an arbitrary state in the sub- 
space spanned by \Tq) and \S). The probability to find 
the electron spins in a state \B) at i > is given by the 
correlation function: 



CsAit) = J2 Pii"") I ("I ® (^1 \A) (37) 



where pi{n) = J^iPiWlil'^ ■ The correlation function has 
the following symmetry: CsAit) = CAB{—t)^ and if \B) 
and \D) are orthogonal states we have CsAit) = 1 — 
CoAit)- Further, we may decompose Cba(0 i^^^o the 
sum of a time-independent term and an interference 
term C)^^(t): 



(38) 



where the overbar is defined in Eq. (|10|l . 

We have further Cg^ = Cba{6K) = Cba{x). Per- 
forming the continuum limit as described in Eq. we 
obtain for the correlation function 



GBA{t) 



dxpi-ao.xoix) {Cba{x) + CgAix^t)) 



C 



BA 



y^int 
^BA 



(t). 



(39) 
(40) 



Here, C^^j is the assymptotic value of the correlator 
CBAit) for t oo. 

We have calculated correlation functions for the fol- 
lowing states: \S) |r"^ = -1) , |To) = +1) , 
\X) ^ |r- - +1) = 71 i\To) + \S)), \Y) ^ Irv = +1) = 
71 (I To) + * I'S'))- The frequency in the interference term 

is always given by s{x) = \/j^ + 4x^. In Table we 
list the integrands according to the notation in Eq. 
(|39|l . From the Heisenberg equation of motion we find 
= — Jt^. which leads to relations for the corrcla- 

dt 

tors. In the notation used in Table we obtain = 
—J [Cyx ~ 5)1 which is satisfied by the results shown in 
Table HJ Similar relations can be derived for the other 
correlators and used to check the results in Tabled We 



CBA{t) 


Cba{x) 




CToSit) 






CTox{t) 


1 1 Jx 

2 s(xy^ 


-IR7Cos(s(a;)t) 


CT,Y{t) 


1 

2 


^ smisix)t) 


Cvxit) 


1 
2 


2iF) sin(s(a;)t) 


CYY{t) 


1 
2 


i COs{s{x)t) 


Cxx{t) 


1 1 2a;2 

2 7{x)^ 


^^cosis{x)t) 



TABLE I: Functions Cba{x) and Ca^(2;,f) according to the 
notation of Eq. 11391 for different correlators (with s{x) = 
V J'^ + 4a;2). Cxx{t) is a linear combination of other correla- 
tors. 



see that Cxx (t) is a linear combination of other correla- 
tors: Cxx{t) = CvYit) + CTos{t). For CtoX and CtoY 
the interference term is an odd function in x. Thus, 
the time dependence vanishes for xq ~ and we have 
CtoX = CtqY = 1/2 for all t. In general, the integral in 
Eq. H39(l is difficult to solve exactly. Thus, we concen- 
trate on several interesting limits. We illustrate this for 
the case of Cyx [t] and give results for the other correla- 
tors. We have 



C'Yxit) = ^+Ini 



^YX 



Cint 
YX — 



Pl;ao,xo (x) 



J 



2s(x) 



s{x)t 



(41) 
(42) 



In the regime of |a;o| 3> ctq the main contribution to the 
integral comes from a narrow region around xq and we 
may approximate ^/(x) ~ where luq = s{xo) and 



too + ^{x - Xo) + 



2s(x) 

in the frequency term s{x) 
For this to be a good approximation, we require (x — 
xoYt ^ 1. We use {x — xq)^ ~ <Tq and thus obtain for 
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the correlator and the range of vahdity in this hmit 



(43) 



LOq 



iUo = ^J-P+^xl (44) 

, , (j2+4a;2)3/2 ^ ^ 

|xo|»ao, t^^^j2^—- (45) 

The results for the other correlators are (with the same 
range of validity) 



0/ cos(tJoi), (46) 



C'^Mt) = -"-4e~'\^oJ cosM, (47) 



Ci?VW = — e"^V*oy sin(c.ot), (48) 



(49) 



CyyW = 1^ Hty cos(cuoi). 



In this limit we obtain a Gaussian decay for all correlators 
on a time scale t'J. = -jr^, — which grows with the absolute 
value of the exchange coupling |J| and with I/ctq- The 
long-time saturation value is 1/2 for Cyx- For some 
of the other correlators we find non-trivial parameter- 
dependent saturation values. In the limit of \xo\ ^ cfq 



we obtain these correlators by the same approximation as 
for the interference term, i.e. we set CsAix) = CsAixo) 
and obtain 



(50) 



^0 

Jxq 



0. 2 ^ l-o|»ao, (51) 



(52) 



For large J the saturation value is quadratic in xq/J for 
CtqS a-nd linear for CtoX- The saturation value for C'tqS 
goes to zero for |J| ^ \xo\ and for CtqX approaches 
1/2. reaches extrema equal to ^ + isign(Ja;o) for 

|J| = 2|xo|. 

Next we consider Eq. for \ J\ max(|a;o|, tJo) and 
find 



s{x) = VJ^ + 4x2 ^ I j| ^ 



2x2 



(53) 



J 



J 



2.(.) 277^ " ^^^^^("^G-^i '- ^''^ 



For Eq. 153|) we have the additional requirement that 
t ^ 7^ — M-1 — lY • Under these approximations we find the 
following result: 



CV"^(t) = sign(J)Qe(i)-5l^'(t)-5|e'(0)exp(^*|J|t-^(l-e'(t))), (55) 

/ t\^^^^ IJI UP 

C(i)^(l-.^j , t^ = i^, |J|»max(|.o|,ao), (56) 



At short times we expand S,'^ (t) ^ 1 + ijr — (^jr^ ■ Keep- 
ing only lowest order in t/tQ in the prefactor and second 
order in the frequency term we obtain 



C^"^(i) = sign(J)ie Hiy sin(c.^i), (57) 



IJI 



, c.^ = |J| + M±^(58) 



t<io = T4' |J| >max(|xo|,cro)- (59) 



4|xo|cro 

4a, 







The |xo| 3> (To limit of this result agrees with the 
I J| ^ |a:o| limit of Eq. (|43|l . Again, we have a Gaus- 
sian decay on the same time scale t'^ as in Eq. H43|) 
{ujo = ^ J2 -f 4a;g ~ | J| for | J| |a;o|)- One interesting 



feature of this correlator is the fact that there is a change 
of phase by tt when the sign of the exchange coupling J 
changes. This feature offers the possibility of measuring 
J even for small values of J through a measurement of 
this correlator. We also list the other correlators in this 
regime: 



^TaX (^) 
^TqY (^) 



J2 



2(x2 + a2)^ H'o/ cos(c.;,t), (60) 



H*o/ cos{Lo'^t), 



—e V J sm(tJot), 

1 'V 
-e 'V^/ cos(wot). 



(61) 
(62) 

(63) 
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Finally, we are also interested in the behavior for large 
t. Thus, we expand Eq. for large times ^{t ^ ip) ^ 
e"^*\/%/t and obtain 



C'^'At) ^ sign(J)e 



J|sin(|J|< 



4(70 1 5 



(64) 



I. /I 



t^t'o^ -r^^ \J\ > inax(|a;o|,cro). (65) 

4(7 n 



For the other correlators we find 



_^cos(|J|i+3^) 



4CTr 



(66) 



C^-^(t) 



-sign(J)e -"0 — ^- ^(67) 



^%a;oV|J|sin(|J|t- 



V|J|cos(|J|t+f) 

1 —■ 



(68) 

(69) 
(70) 



Thus, the transverse components of the pseudo-spin have 
a slower decay (~ t~^/^) than the longitudinal compo- 
nent t^'^/'^). This results from the fact that the Hamil- 
tonian only has fluctuations along only one direction. 



V. ANALYSIS OF VSWAP 



In this section we analyze the VSWAP gate using 
the correlation functions d erived in the previous sec- 
tion, i.e., we analyze the vSWAP gate taking into ac- 
count the hyperfine-induced dccohcrcnce. The VSWAP 
gate and single-qubit operations can be used to perform 
the quantum XOR gate (CNOT) which, in combination 
with single-qubit operations, is sufficient for universal 
quantum computationiil In Ref. |^ implementation of 
vSWAP has been demonstrated. However, in these ex- 
periments there was a contrast reduction of ~ 40%. Here 
we show that taking into account hyperfine induced de- 
coherence, still near-unit fidelity can be obtained for this 
operation. 

The Hamiltonian of Eq. |^ induces unitary time evo- 
lution on the states of the system: |'0(i)) = U{t) |V'(0)) 
with U{t) = Tcxp(-i ^lH{t')dt'). We assume that J 
and xq can be switched adiabatically^ on a time scale 
that is much shorter than the time required for the gate 
operation and thus the time evolution operator at time 
Tc has the form 



Us = exp {-iTsH) 



(71) 



In a Bloch-sphere picture this operator induces a rota- 
tion about an axis in the plane spanned by eigenstates 
of T- and r^ \X) = and \S) = (|n) - |iT))/%/2^ 



The axis of rotation is determined by the parameters J 
and ccq. Through such an operation any state may be 
rotated into any other state on the Bloch sphere. Thus, 
it is possible to rotate from \S) to any initial state in 
the subspace of S*^ = by a single operation. This is 
important since initialization to the singlet is feasible by 
preparing a ground-state singlet with both electrons on 
the same dot and then changing the biasnS. We now inves- 
tigate initialization to the state \X) taking into account 
hyperfine-induced decoherence. The scheme we propose 
here is different from the one used in Ref. |^, where adi- 
abatic passage from the singlet to the |ti)-state is used. 
Our scheme requires control of xo. We assume the system 
to be in the singlet state \S) at t = and then switch J 
and xq such that J — — 2x0 and |a;o| 3> ctq. In a Bloch- 
spherc picture, this corresponds to a rotation about an 
axis that halves the angle between \S) and \X). Since 
Cxs{t) = Csx{-t) = 1 - CTox{-t) we have, for the 
above choice of parameters, according to Eqs. l|T7|l and 



Cxsit) 




This correlator reaches its maximum for A/2|t/|i = tt, i.e., 
at Ts = . The time scale for the Gaussian decay is 

t" — — . To approach unit fidelity we therefore require 
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I J| ^ iJoj which is the case in the range of validity of the 
above correlator since |xo| 3> ctq and J and xq are of the 
same order. At t = we switch J to zero and since 
\X) ®\n) is an eigenstate of the remaining Hamiltonian, 
the system remains in this product state, untouched by 
decoherence induced via the nuclear spins. This scheme 
thus provides a way to initialize the double quantum dot 
system to the state \X) = ^(iTb) + \S)) = |tJ,), where 
arrows denote the z-component of the electron spin in 
each dot. In the same way, it is also possible to initialize 
in the state \~X) = |r- = -1) = ^(|To) - \S)) = 
by switching to J = 2xq. 

It was already proposed in Ref. to implement the 
VSWAP gate by pulsing the exchange interaction J be- 
tween the two dots. Here we give a detailed analysis of 
the VSWAP gate taking into account hyperfine-induced 
decoherence. 

The SWAP operation acts on the basis of the two- 
electron system as: ^ Ui> , UT> ^ ITi),ITi) ^ 

lit) , ITT) |TT>- The SWAP is an operation that acts 
only on the subspace of S*^ = and leaves the states 
Itt) and unchanged. In the system we consider 
this is naturally implemented through the large Zee- 
man splitting that separates |tt) and \[[) from the sin- 
glet and the S"^ ~ triplet. In order to analyze the 
SWAP in the 5*^ = subspace wc consider the regime 
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of \J\ » max(a;o, (To)- The correlator C-x,xit) gives 
the probabihty of bemg in the state \—X) = jjl) for 
a system initiahzed in \X) = Due to the sym- 

metry relations for the correlation functions we have 
C_x,x(i) = 1 - Cxx{t) = 1 - CYvit) - CTosit) and 
thus find (using Eqs. (|Hn|l and and neglecting terms 
of order (ctq + x'^)/ J'^), 

C.x,x{t) = 1 - Cxxit) « i - ^e"Kt) cos(| J|0, 

(75) 

t'o = *f | , |J| > max(|a;o|,fTo), i<t^ = -|4- 
4(70 Fo I 4crQ 

(76) 

Wc obtain the maximum value for this correlator when 
r, = rrr • The Gaussian has a decay time of t'J. = . , , 

^ \J\ J V 4(To|2;o|' 

SO for xq ^ the Gaussian decay is negligible and we 
obtain unit fidelity for this SWAP operation |tj) |jt) 
up to a global phase factor (which is not visible in the 
correlator). 

Prom the SWAP operation it is only a small step to- 
wards the vSWAP which we obtain when we let the sys- 
tem evolve with the same parameter values but for only 
half the time. Starting in the state \X) we obtain \Y) 
after applying a VSWAP. Por large | J| we find for the 
correlator Cyx in the limit xq ^ 

Cr;f(t) = i+sign(J)ie^Kt) sin(|J|t), (77) 

*o = Z^T^' l-^l »max(|xo|,tTo), i<to = ]4- 

(78) 

Here again the time scale of the Gaussian decay is 
and approaches infinity for — *■ 0. The time during 
which we have to operate with these values of the param- 
eters J and xo is now Ts = 2[7\- calculations show 
that for the time during which J is pulsed high there is 
a regime in which unit fidelity may be approached. The 
reduced visibility in the experiment^ may be due to sev- 
eral reasons such as reduced visibility in the readout of 
lit) or the initialization of |tj). 

VI. CONCLUSION 

Wc have developed a method that uses the measure- 
ment of a Rabi resonance in the quantum-dot spin qubit 
to narrow the distribution of the nuclear spin states. This 
method relics on Rabi oscillations induced via an oscil- 
lation of the singlet-triplet splitting J in the subspace 
5"^ = of two electrons in a double quantum dot forming 
a two-qubit system. Purther, we have calculated several 
correlators in the S*^ = subspace for static J and found 
that the transverse components of pseudo-spin have a 



slower decay than the longitudinal one. We have also dis - 
cussed the implementation and fidelity of the VSWAP- 
gate in this system and the initialization to the jtj), |],|) 
states. 
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APPENDIX A: DRIFT IN 

In addition to spin diffusion, driven by the nuclear 
dipole-dipole interaction, there may also be a change in 
Sh^ due to corrections to the projected effective Hamil- 
tonian considered here (see Ref. Iisl Appendix B for 
details). After tracing out the electron pseudo-spin in 
state ps, these correction terms give rise to an electron- 
mediated nuclear spin-spin interaction which, in general, 
takes the form of an anisotropic (XYZ) Heisenberg inter- 
action 

Here, the indices i and j run over all nuclear spin sites. 

We use the corrections to leading order in the inverse 
Zeeman splitting 1/ez (e^ = q^bB) given in Ref. IT^ 
This gives the typical value of the exchange constants 
~ A^/N'^ez- Assuming an unpolarized nuclear spin 
state, each nuclear spin will therefore precess in an effec- 
tive mean field generated by all other spins in the dot of 
typical magnitude 

/leff- Va^IJ^I ^ AViVie,. (A2) 

This effective field will result in precession of the nuclear 
spins about an arbitrary angle (and hence, may change 
the value of (5/i^) on a time scale 

Tp - N^e^/A^ ~ 10-2 s, (A3) 

where we have assumed N = 10® nuclear spins within 
the quantum dot. and ez/gpB = A/gp,B — 3.5 T for the 
time estimate. This is only a worst-case estimate, which 
neglects the effects, e.g., of a Knight-shift gradient (due 
to strong confinement of the electron) , which may further 
weaken the dynamical effect discussed here. We expect 
the dipolar nuclear spin diffusion time to be the limiting 
time scale for nuclear spin dynamics, in light of exper- 
iments on diffusion near donor impurities in GaAs^ If 
the effect giving rise to Tp in Eq. ljA3|l were significant, it 
could be further suppressed by choosing a larger quantum 



12 



dot size or stronger magnetic field, thus allowing many 
electron spin measurements on the time scale of variation 
of 6h\ 



APPENDIX B: MEASUREMENT 

In this appendix we describe how a single measure- 
ment of the two-electron system affects the nuclear spin 
state. We give the analytical expression for the diago- 
nal elements of the nuclear spin density operator after a 
measurement. 

At t=0 the system is described by the following density 
operator 



with 



p(0)=Pe(0)®p/(0)= 1+) 



.^p,|V^l)(^l|, (Bl) 



with nuclear spin state = X^n '^n !"■)■ The Hamil- 
tonian Hq of Eq.(^ acts on the the nuclear-spin sys- 
tem as Ho \n) — Hn \'n), where in i/„ the operator 
has been replaced by 5h^ (because 5h^ \n) = 5h^ |n)). 
Since [HQ,5h'] — 0, only the diagonal elements of the 
nuclear density operator pi (in the basis of 6h^) enter 
in matrix elements for operators acting only on the two- 
electron system. As described in Section FlI Al these di- 
agonal elements pi{n) = pi{n, 0) = (7i| Tre{p{0)} \n) de- 
scribe a continuous Gaussian distribution in the contin- 
uum limit. The trace over the electron system is defined 
as Trepit) = (+|p(i)|+) + {-\p{t)\-) and for p/(n,0) 
we have 



pi{n,0) = 



,i |2 



(B2) 



The time evolution operators U{t) and J7„(i) are defined 
through iij{t) = Ho{t)U{t) and iUnit) = iJ„(i)f/„(i) and 
thus the density operator p{0) evolves under the Hamil- 
tonian Hq as 

p{t) = uit)pio)uHt) 

= U{t) Le(0)®^^Ka>r|n)(?| j UHt) 

Y i n,l j 

= [f/n(t)Pe(0)C//(t) ® ^P.a>r \n) (Z| J . 

(B3) 

At time tm a measurement in the basis of |-|-) and |— ) is 
performed on one single two-electron system coupled to 
nuclear spins. Since the outcome of this measurement is 
known, the state of the system after the measurement ia^ 
(the result depends on whether |+) or |— ) was measured) 

|±) (±|p(UI±)(±l 



(tm) 



^ (|±) (±1 C/„(i™)pe(0)f//(t„) |±) (±1 



)^Ka>f \n) (;|) — i 



P±(i™)' 



P±(t,„) - TriTr,{\±){±\p{Ua)} 



(B5) 



where TrjA = X]n("'l^l"') ^^"^ is given in 

Eq.Q. Here, P^{tm) is the probability to measure |±) 
at time tm- We are mainly interested in the diagonal el- 
ements of the nuclear density operator pj after the mea- 
surement. 



p[^'±)(n,i„) = {n\Tr,p^''^\tra)\n) 

= |^|^(±|[/„(t„)p,(0)[/t(t„)|±) 



P/(»,0) 1 



(l±(T^(i„0)„) 



(B6) 



Using Eq.Q we find 



p\^''^\n,trn) 



Pi{n,0) 1 / 2(a, - Lu'f 
P+(t,„)2 [{n„ -Lor + {j/2)^ 
{j/2)^{l + cosiuj't^)) 



(B7) 



and 



(i,_) piin,0) l {j/2y{l~cos{cj't^)) 



where oj' is given in Eq. jSJ and depends on the eigenvalue 
of the nuclear spin eigenstate through ri„. 

Parenthetically, we note that in the case (not described 
in this article) where the measurement is performed on 
an ensemble of many different double quantum dots, the 
state of the ensemble after the measurement is^ 

+ \-){-\Un{tm)PemU}{Ua)\+) (+| 



(B9) 



and the nuclear-spin distribution has not changed. If a 
complete measurement of the Rabi-resonancc lineshape 
would be performed on an ensemble of double dots, the 
result would be the Voigt profile described in Sec. Ill Al 



(B4) 
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